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Abstract 

The motion of single quantum particle through Bose-Einstein condensate (BEC) is considered 
within perturbation theory with respect to the particle-BEC interaction. The Hamiltonian of BEC 
is diagonalized by means of Bogoliubov's method. The process of dissipation due to the creation 
of excitation in BEC is analyzed and the dissipation rate is calculated in the lowest order of 
perturbation theory. The Landau's criterion for energy dissipation in BEC is then recovered. The 
energy spectrum of the impurity particle due to the interaction with BEC and its effective mass 
are evaluated. 
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I. INTRODUCTION 



Recent experiments on the rotational motion of molecules in the superfluid helium 
droplets gave rise to a new interest in superfluidity on a microscopic scale p, |7| . A unique 
properties of the measured rotational spectra pose a large number questions, such as: what 
is the collective molecule/superfluid wavefunction that describes sharp rotational states ob- 
served in experiments; what are the properties of finite systems and how is the limit of a 

microscopic phenomena were studied in the BEC 



of Sodium atoms in magnetic traps 



231 ] . A linear motion of impurities was shown to be 



dissipationless for the speeds below the condensate speed of sound. A large number of 
theoretical works described the molecule/ droplet s yste m using imaginary time path inte- 



gral Monte- Carlo approaches 



0,Q. 
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12, m 



While in certain cases remarkable 



151 ]. those works are strictly limited 



agreement with experimental constants was obtained 
to the calculation of statistical properties and thus provides no real understanding of the 
microscopic nature of the dissipationless motion. The latter can only be established by con- 
sidering a real time dynamics. Several theoretical works considered a motion of impurity 
through dilute BEC. Recent works considered a macroscopic particle interacting with dilute 
BEC [3j. In this case the motion of particle is equivalent to the BEC in a time dependent 
external potential. This problem was treated by solving time-dependent Gross-Pitaevskii 



equations 



4|. A microscopic particle interacting with the dilute BEC in Bogoliubov's ap- 



proximation has been considered by several authors using general Golden rule considerations 



19l . l20l . l2lj . These works were based on the result of Miller et al 



22 1 which was obtained 



using time-independent perturbation theory. The Bogoliubov's treatment has also been suc- 
cessfully used for the investigation of the force acting on the impurity particle due to the 
quantum fluctuations in BEC 



24 



251 ]. Some authors treated a partic 



16 



e str ong ly interacting 



13, 



with Bogoliubov's BEC and found a possibility of self localization 

In this work we attempt to describe impurity moving through BEC as a microscopic 
particle within time-dependent perturbation theory. We introduce the Hamiltonian of a 
quantum particle moving within the interacting Bose gas. No assumption is made about the 
relative mass of an impurity compared to that of the Bose particles. We present the analytic 
results for this problem in the limit of the dilute Bose gas at zero temperature. After the 
introduction of the general Hamiltonian, the Bogoliubov's approximation is made to convert 



2 



the Hamiltonian to the diagonal form. Then the problem is reduced to the quantum particle 
moving through the gas of non-interacting Bogoliubov's excitations. The initial conditions 
are chosen such that an impurity is not correlated with the BEC. The static properties are 
achieved after the system reaches thermal equilibrium in a long time limit of our treatment. 

We obtain the expansion for the transition amplitude in powers of particle-BEC inter- 
action and restrict ourselves to the calculation of the leading term of this expansion that 
corresponds to the Golden rule approximation. The latter allows us to calculate various 
observables such as: state-to-state transition probability, dissipation rate, and the effective 
mass of the particle. Then we obtain the criterion for the critical momentum of a particle 
below which its motion becomes non-dissipative, i.e. the microscopic version of the Landau's 
criterion. This criterion comes from the exact momentum-energy conservation law. Finally, 
we calculate the branch of the excitation spectrum corresponding to the moving impurity, 
i.e. the energy of impurity as a function of its momentum. Expanding the energy at small 
momenta it is found that the linear part of energy vanishes and the impurity behaves like 
a free particle with an effective mass that is determined by the second order derivative of 
energy with respect to momentum. We show that our results agree with a time-independent 
treatment of Miller within long time limit. 

II. MODEL HAMILTONIAN 

In this section the model for the interacting Bose gas is discussed as well as the approxi- 
mations that are made in our treatment. The Hamiltonian of interacting Bose particles in 
coordinate representation is given by 



Here and throughout the paper we set h — 1. In the secondary quantization the above 
Hamiltonian reads 




(1) 




(2) 



p 



P1.P2.P 



where U (p) is the Fourier transform of the interaction potential 




(3) 
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We will concentrate on the case of dilute gas r <C ri -1 / 3 where r is the range of potential 
on which U (r) differs from zero significantly and n denotes density of gas. So the integral 
([3]) is significant in the domain p < 1/tq while characteristic momenta of particles are of 
the order of p c < n 1 ^ 3 . Thus due to diluteness condition p p c , and the Fourier transform 
of the interaction U(p) in Eq. (j2j) can be replaced by Uq — J U(r)dr. The zero Fourier 
component Uo is connected with length of s-scattering in first order Born approximation as 
following 

a, = £f/ . (4) 

Let us consider degenerate Bose gas at zero temperature. In this case the Hamiltonian (J2J) 
can be reduced to the diagonal form with the help of the Bogoliubov's method Q]. The main 
assumption of the original treatment is that most of the particles of the degenerate gas stay 
in the ground state with zero energy and the number of particles with non-zero momentum 
iVp^o is much less then Nq = N p= q. Thus the creation/annihilation operators b p=0 and 6 p= o 
can be replaced by the macroscopic number y/No. So one can expand the interaction part of 



the Hamiltonian (J2J) in powers of y/N p ^o/N an d leave the terms of expansion up to second 
order only. After the described procedure the Hamiltonian of the Bose gas attains the form 



p 



iV 2 - N + N (K b - P + M-p + 2 K b r>) 
p^o 



(5) 



where A" = N + N pj i Q . The next step is to diagonalize the Hamiltonian (jSJ) introducing new 
bosonic creation/annihilation operators B£ and B p related to old operators bi~ and b p by 
linear transformation 

b P = a p B p + /3 p B+ p , b+ = a p B+ + (3 P B_ P . (6) 

Substituting Eqs. ((BJ) into the Hamiltonian (jSJ) and requiring the Hamiltonian to have diag- 
onal form, one gets for the transformation coefficients |l| 

Qp = / 2 1 ' = / 2 1 ' 

VMp - 1 V^p - 1 
e(p) +p 2 /2m + nU 
= Tj ' 

together with the Hamiltonian in new representation 

H B = E +Y,<P)B;B P . (8) 
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The excitation spectrum e(p) of the system yields 



p2 / p2 



and has the phonon-like behavior at low momenta, i.e. e(p — > 0) = p^nU^jm = pc, where 
c is the speed of sound. The ground state energy of BEC is given by 

Now let us introduce a single quantum particle with mass M and coordinate R interacting 
with the environment of Bose gas discussed above. The whole system is then described by 
following Hamiltonian 

h = -w,^ + Hb + h - <"> 

with the particle-environment interaction 

H * = v T, b U b * e ~ ikR - ( 12 ) 

p,k 

Here the coupling constant g is determined as zero Fourier component of the system- 
environment interaction. In the spirit of Bogoliubov's theory we expand the interaction 



Hamiltonian (fl2l) in powers of yiVp^j/iVo and leave the terms of zeroth and first order, i.e. 
leading terms. Then the Hamiltonian ( fl2l) takes the form 



Hi = gn+ 9 -^J2 ( & p e " ipR + Ve ipR ) • (13) 

We neglected the terms of the order of N p ^q in (fT5|) . As it becomes clear below, the con- 
tribution of the higher order terms to the dissipation process is of the order of N p ^q which 
would be an excess precision compared with the Bogoliubov's approximation made above. 
Expressing operators 6+ and b p through new operators B£ and B p , and using the secondary 
quantization representation for the system particle, for the full Hamiltonian ( TTT1) one gets 

q 2 

H = K + Yl 2M a q flq + 6 (p) B p B p + 7 p ( a J- P a q fl P + a q+ P a <i B p) > ( 14 ) 
q p q.p^o 



g / Np 2 
7 P ' ' 



V y 2me(p) 

Here operators a+/a q create/annihilate the particle in state |q), and the ground state energy 
E' Q = Eq + gn is shifted with respect to the E due to the system-condensate interaction. 



III. TRANSITION PROBABILITY AND DISSIPATION RATE 



The next step is the calculation of the transition probability for the whole system from 
some initial state \i) to final state |/). So our aim is to construct the perturbative expansion 
for the matrix element of the evolution operator a,/(t) = (f\e~ lHt \i) in powers of particle- 
BEC interaction (fl3l) . Here we restrict ourselves to the calculation of the leading, i.e. lowest 
order term of this expansion only. Let us write the matrix element of evolution operator in 
interaction representation 

<Hf(t) = (f\e- lHot Texp f-i J F/(s)ds^ \i) , (15) 
H = H-E I , ifj(s) = e iHoS H ie - lHoS . (16) 

Initial and final states of the system particle are supposed to be free particle states |q) = 
_^ e «qR with momenta and q/, respectively. The Bose environment is initially in its 
ground state |0)b and then evolutes into some final state \f)s- Note that the state \0)b is 
the ground state of the diagonalized Hamiltonian (jSJ) with the ground state energy E' and 
it does not equal the state in which all of the Bose particles have zero momentum. Below 
we will find the transition amplitude as an expansion in powers of coupling constant g, i.e. 

Oif(t) = af/(t) + + + - , a$(t) ~ g n . (17) 

The zeroth order term is given by a$(t) = 5 q/>qi exp [-i(E' - qf/2M)t). The first order term 
of expansion describes the process of creation of excitation in BEC with the momentum pf 
due to the interaction with the particle 

ft 

i \ .. /„ I /. i -illnU-."\ + r>+ -ill„.- \ .. \ li'i^ 



og^t) i / ^ETp(P/lB(q/|e- Wo(t - s) aJ_ p a q Sj e -^|q i )|0) J 

J ° P,q 

= -ilpA m i-Pf I dsex P [-i(E -q'>/2M)t-iL J (p)(t-s)] , (18) 

where w(p) = e(p) +p 2 /2M — qjp/M . The term of second order is defined by two processes: 
first, is creation of two excitation in BEC; and second, is creation and annihilation of one 
excitation in BEC while the particle finally remains in its initial state. We will consider the 
amplitude of the latter process since only this amplitude gives the contribution of second 
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order into the transition probability. Thus the second order matrix element reads 

ft rt 



\ [ ds [ d S /(0|B(q/|e-^T(F / (5)^/(^))|q,)|0)i ? 
* Jo Jo 



= -SnfM I ds I ds ' Q ( s ~ s "> exp H ( E 'o ~ <li/ 2M ) 1 ~ iuJ (p)( s - s ')] ■ ( 19 ) 
p Jo Jo 

Finally, one can find the probability of the particle transition from state |q$) to |q/) in the 
following form 



<u,if 



7 C 



dse~ isu}(p) 



I 2 



(20) 



The standard way to deal with the integral J * dse~ luJS in the expression for the transition 
probability fTSOj) is to evaluate its value in the limit of long time interval t ^> 1/u. So for 
the square of the absolute value of this integral one gets 

2 

— ► 2nt5{u(p)) . (21) 



dse~ lsuj{p) 



t— »CXD 



First, let us discuss the probability of the particle of leaving its initial state |qj) and 
moving to the state with a different arbitrary momentum |q/) in a unit of time with creation 
of the Bogoliubov's excitation in the surrounding BEC with the momentum p = q/ — q^- In 
the long time approximation the probability reads 



dw(p) = 2njl5(u(p))dt 



(22) 



So the maximum momentum pm that can be transmitted from particle to BEC is collinear 
to qj and determined by the crossing point of the line qip/M and the curve e(p) +p 2 /2M 

2 



Pm 



y/q 2 c + (M/m)*(<g-q, 



c 2 ),- 



q c = Mc. 



(23) 



1 - (M/m) 2 

Besides, the maximum angle 8m between vectors q« and p is given by the condition cos 6m = 
q c /qi- Thus we have the Cherenkov like emission of the BEC excitation with the length of 
the momentum < |p| < pm into the cone with the angle 9 m along the vector q^. 

Next, we can define the transition rate as the probability of the particle to leave its initial 
state |qj) in a unit of time as 



E 



dw(u(p)) 



(24) 
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In the thermodynamic limit (N — ► oo, V — > oo, iV/V^ = n) one has to replace the sum 
over momenta by the integral / <^P i n t ne right hand side of Eq. and for the 
transition rate one gets 



2 p 1 

ng p 



Fr = / dp-?y MP)) • (25) 

S^m J e[p) 



Let us write the integral in Eq. ( |25|) in the following form 



poo p 3 p1 

/ dp—— / cfe<5(x — xo) , (26) 
Jo e iP) Vo 



47rmgj Jo 

where x = pqi/TO and the value Xq is determined from the condition uj(p) = 0. So the 
integral over x is unit only if xq < 1 , or 

e(p) + p 2 /2M < pqi/M . (27) 

Clearly the condition (1271) is never fulfilled if the absolute value of the initial momentum of 
the particle qi is less then Me(p — > 0)/p = Mc = p c , and then the integral in the expression 
for the transition rate IV is always zero. At this point we recovered the Landau's criterion 
- there is no energy dissipation if the speed of particle = qi/M moving through the 
interacting Bose gas is less than speed of sound. 

To evaluate the rate for momenta above the critical one we proceed further with the 
integral ( 1261) 

nMg 2 [*" p 3 
T T = / dp— - . 28 

Here the integration is performed over momentum p from zero to the crossing point p = Pm 
given by Eq. (|23|) . Finally, the transition rate reads 

Mmng 2 



litqi 



( \ 2 1 / i i e(p M ) +P M / 2m 
e[PM) - mc In ( H — 



(29) 



Let us consider the limit case when the initial momentum q$ is close to its critical value q c . 
Expanding the expression (J23l) in powers of small — q c and assuming the value Mj m to 
be finite, for the upper limit of integration we obtain pm ~ 2(gj — q c ). Then expanding the 
right hand side of Eq. ( |29l) in powers of pu / mc, for the leading term of expansion one gets 



2na 2 

T T {q l -qc) = ^ JL - 2 {q l -q c f ■ (30) 
SiTmc z 



So the lowest term of expansion appears to be of the third order of small value — q c . In 
the opposite case of large initial momentum ^> q c one can obtain the following expression 
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rx = ^f^| , (3D 



for the transition rate 

= V*L 

which is surely independent on speed of sound. 

Our next example is the calculation of the energy transferred from particle to BEC as a 
function of time 

E(t) = 5>(P)7S fdse—^ = 2ntjf - [ dp Tp e(p)5(u,(p)) . (32) 

Here we can determine the dissipation rate as the amount of energy transferred in unit of 
time 

r E = m. = -g- I dp ^M P )) . (33) 

dt oTt 171 J 

Calculating the integral in ( 1331) in the same manner as described above, for the dissipation 
rate one gets 

T E = ¥%&L . (34) 

In the limit case of massive particle M ^> m the expression fl34|) converges to the following 
value 

Fi? = ^ • (35) 

This result exactly corresponds to the rate of energy dissipation obtained in [3J for the 
classical particle moving through BEC with constant velocity. 



IV. PARTICLE ENERGY SPECTRUM AND EFFECTIVE MASS 

In this section we will consider the non-dissipative motion of the impurity though the BEC 
when the initial momentum of the impurity particle is less than the critical momentum 
q c = Mc at which the dissipation process occurs according to the Landau's criterium. In 
the absence of dissipation one would expect that the impurity will move as a free particle 
but its energy should be shifted due to the coupling with the BEC environment. Thus our 
purpose is to calculate the contribution to the energy of the impurity due to the interaction 
with BEC and the effective mass of the particle. 

Let us start with the calculation of the full energy of the particle-BEC system as a 
function of time. Following the perurbative treatment developed above, one can write the 
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energy of the system as an expansion in powers of coupling parameter g 

E(t) = E ( -°\t) + E (1 \t) + E {2) (t) , 



(36) 



coming from the expansion for the evolution of density matrix 

p(t)=pW(t)+p^(t)+pM(t). (37) 

The system energy as a function of time is then given by E(t) = Tr (Hp(t)), where H is the 
Hamiltonian of the full system (fT4l) and Tr denotes the trace over all possible states. The 
zeroth order contribution corresponding to the unperturbed system is 

= Tr (H oP ^) =E' + ^. (38) 

The first order does not bring any contribution since Tr (H p^ (t) + Hjp^(t)^j = 0. The 
second order term consists of two terms E^ = Tr (H p^) + Tr (Hjp^). Let us consider 
the first term 

IV (H 0f P) = £ (if + <P,) + K) + <U (ii + K), (39) 

where the transition probability w\J is defined by Eq. (1201) as w\j = Wif — 1. We are 
interested again in the long time limit when the whole system reaches its equilibrium state, 
so the upper limit of integration over time in the right hand side of Eq. ( 1201) must be replaced 
by infinity. Since here we consider the case of dissipationless motion of the impurity, u>(p) ^ 
and instead of the limit (1211 one has to use following formula 

2 



t 

dse- isw{p) 



1 



(40) 



Hence for the term ( |39|) we have 

W^E^y (41) 

p 

The remaining term Tr (Hip^*) is given by 

Tr (H IP U) = J2 (</|F/|^> (1) (*)l/> + cc) . (42) 
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Using the expression for the interaction Hamiltonian Eq. f[T3"j) and noting that (i\p^(t)\f) = 
af) (t)aif (t), for the above term one gets 

Tr (H IP U) = -2 V 7 p [ ds sin(^(p)) . (43) 
p Jo 

In the long time limit J °° ds sin (su(p)) = V^^, where V denotes principal value. Finally, 
the energy shift given by the sum of terms fl39l) and (j33J) reads 

E {2) = -T4^ = ~^K- f dp , f, , • (44) 



p 



As was mentioned above, here we are interested in evaluation of the energy of the impurity 
particle in case of relatively small initial momentum qi. In this case the energy can be 
written as an expansion in powers of small parameter qi/q c up to second order, i.e. 

E p ( gi ) = E p (0) + ^- . (45) 

It turns out that the term linear in disappears, and the effective mass of the particle in 
(1451) is determined as follows 



1 1 d 2 E^ 
+ 



M ef M ' dqf 



(46) 

qi=0 



Let us start with the calculation of the zeroth order contribution in powers of qi 

E(0)=n - — ^ r P ' dp (47) 

pU ng 4tt 2 m J e(p)(e(p)+p 2 /2M) ' 1 ' 

The integral in the right hand side of Eq. f|47|) diverges at large momentum like 4mm r J °° dp, 

where m r = (l/m+ l/M) -1 is reduced mass of the system of the relevant and the BEC 

particles. The standard way to avoid the such kind of divergence is to renormalize the 

coupling constant introducing some another expansion parameter. Let us determine the 

coupling parameter g as an expansion in powers of scattering length a 

2vra . . 

g = (1 + kcl + ...) , (48) 

m r 

and then remove the divergence from the expansion for the energy ( T4T1) to the expansion 
PU) . Substituting equation (J151) into Eq. (jlT]) . leaving terms up to second order in powers of 
scattering length a and requiring the energy E p (0) to be finite, for the renormalized coupling 
constant one gets 

2na ( 2a f 00 \ 

" = -^{ 1 + tIv- (49) 

11 



The actual physical reason of the divergence is that the zero Fourier component of the 
interaction is used as the coupling constant g. The procedure of renormalization of coupling 
constant corresponds to taking into account the second order Born approximation and can 
also be done employing direct expansion of the scattering length in powers of coupling 
constant Using the relation (j49j) . the energy E p (0) is now finite and reads 



E p (0) 
h{x) 



2ima 



Aarac r ( m 

1 + — 7 » (m 



(50) 



x-\/x 2 — 1 — ln(x + a/x 2 — P 



(x - l)Vx 2 - 1 

Evaluating the term of the second order expansion in the expression for the energy 



d 2 E^ 



dqt 



na 



Qi=0 



p 6 dp 



3 mm 2 r M 2 J e ( p ) ( e ( p ) + p 2 / 2 M) 



3 ' 



(51) 



for the effective mass of the particle we have 

32 



M, 



ef 



M 



na 



3 Mc(m r /m) 



m 
M 



n -1 



where 



h(y) 



;i + 2y 2 ) \n{y + ^y 2 - 1) - 3y^y 2 - 1 



(52) 



(53) 



2(y 2 - l) 5 /2 

The behavior of the functions Iq and I\ in expressions for the energy fl50l) and the effective 
mass (|52|) . respectively, is shown in Fig. ([1]). In case of heavy particle M m its effective 
mass yields M e f = M(l — 8irna 2 /3Mc) _1 , the same result as was obtained in [ij]. 



V. CONCLUSION 



The goal of this paper is the understanding of the effect of superfluidity reflected on a 
motion of a microscopic impurity inside a BEC environment. Quantitatively our system is 
very different from the realistic liquid helium. However, for dilute BEC we obtained physical 
behavior that gives qualitative understanding of manifestation superfluidity on a microscopic 
scale. We calculated the transition amplitude to the second order of the perturbation theory. 
The inelastic contribution to this amplitude responsible for the dissipation comes from the 
first order of the perturbation expansion. It is shown that at the velocities less than speed 
of sound this term is zero and we recover the exact Landau's criterion. The elastic part 
corresponds to the second order term and is responsible for the shift of the particle energy. 
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Divergence of this energy shift which appears because of the delta- functional potential is 
avoided by renormalization of coupling constant. The energy of the particle is calculated as 
a function of its momentum. At small momenta this function is quadratic, thus the impurity 
moves like a free particle but its mass is shifted due to the interaction with the BEC. The 
natural extension of this work is the calculation of the Green's function of an impurity and 
the development of the perturbation theory for its poles. This work is currently in progress. 
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